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Gravitational lensing distorts the cosmic microwave background (CMB) anisotropics and imprints 
a characteristic pattern onto it. The distortions depend on the projected matter density between 
today and redshift z ~ 1100. In this paper we develop a method for a direct reconstruction of the 
projected matter density from the CMB anisotropics. This reconstruction is obtained by averaging 
over quadratic combinations of the derivatives of CMB field. We test the method using simulations 
and show that it can successfully recover projected density profile of a cluster of galaxies if there 
are measurable anisotropies on scales smaller than the characteristic cluster size. In the absence of 
sufficient small scale power the reconstructed maps have low signal to noise on individual structures, 
but can give a positive detection of the power spectrum or when cross correlated with other maps of 
large scale structure. We develop an analytic method to reconstruct the power spectrum including 
the effects of noise and beam smoothing. Tests with Monte Carlo simulations show that we can 
recover the input power spectrum both on large and small scales, provided that we use maps with 
sufficiently low noise and high angular resolution. 

PACS numbers: 98.80.Es,95.85.Bh,98.35.Ce,98.70.Vc 



I. INTRODUCTION 



Cosmic microwave background (CMB) anisotropies have the promise to revolutionize the field of cosmology in the 
following decade. Using ground based, balloon and space experiments we will be able to map the microwave sky over 
a large range of angular scales and frequencies. The expected characteristic pattern of acoustic oscillations generated 
by the primary CMB anisotropies around z ~ 1100 will provide a wealth of information that should constrain many 
of cosmological parameters to a high accuracy Other sources closer to us also contribute to the microwave 
sky. Some of these are the foreground emission from our own galaxy and from the galaxies along the line of sight, 
Sunyaev-Zeldovich emission from clusters, signatures of patchy reionization, etc. 

Another effect that modifies the CMB sky is gravitational lensing. Dark matter distributed along the line of sight 
between z ~ 1100 and present deflects the light and induces distortions in the pattern of the CMB anisotropies. Its 
effect on the CMB power spectrum has been thoroughly investigated ||J| . The conclusion from these works is that the 
lensing effect is a small but non-negligible. On large and intermediate scales lensing smoothes the acoustic oscillations 
, while on very small scales it creates additional power. Both these effects could help break some of the parameter 
degeneracies in the CMB Q| . Although in models normalized to cluster abundances the effect of gravitational lensing 
is small it could be important for Planck satellite mission f|. The lensing effect on the CMB power spectrum is now 
included in the standard CMB primary anisotropies routine ||. 

Gravitational lensing is directly sensitive to the matter distribution up to z ~ 1100, so a detection of this effect 
would provide important information about matter distribution on large scales and high redshifts, which is not directly 
attainable by any other means. Such information would not require any additional assumptions such as how light 
traces mass or how nonlinear structures form. Because of this it is worth investigating if there are other signatures 
imprinted on the CMB which would be more easily accessible to the future observations than the effect on the 
CMB power spectrum. Some of these, such as four-point function and ellipticity distribution of peaks have already 
been explored by Bernardeau Q. These particular signatures of the lensing effect were found to be rather weak, but 
nevertheless marginally detectable with Planck. They could provide additional constraints on the amplitude of matter 
power spectrum. Similar conclusion has also been reached by Kaiser 

The purpose of this paper is to present a new approach to identify gravitational lensing in CMB. The power of the 
method developed here is that, unlike previous attempts, it allows a full 2-d reconstruction of projected matter density 
between us and the last scattering surface at z ~ 1100. Such a map can be used to search for clusters at high redshifts 
or study density fluctuations on the largest scales that are not directly accessible otherwise. It can also be correlated 
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with other maps of interest, such as those of the Sunyaev-Zeldovich effect, X-ray background, weak lensing, galaxy 
clustering and the CMB itself. In this paper we develop the method and test it on simulated maps, showing that 
it can give an unbiased estimate of the underlying projected density field and its power spectrum. Whether or not 
the reconstruction from CMB can be successfully applied to the real data depends on the level of CMB anisotropics, 
angular resolution and noise characteristics of particular experiments as well as the amplitude of matter fluctuations. 
In this paper we discuss these issues in detail and show the conditions that need to be satisfied for the method to 
work in practice. 



II. RECONSTRUCTING THE PROJECTED MASS DENSITY: FORMALISM 



The large scale density fluctuations in the universe induce random deflections in the direction of the CMB photons 
as they propagate from the last scattering surface to us. This effect not only alters the power spectrum of both 
the temperature and polarization anisotropics but also introduces non Gaussian distortions in the maps. The 
quantity responsible for the deflections is the projected mass density or convergence k, defined more precisely below. 
In this section we develop the formalism to measure convergence k based on its lensing effect on the CMB maps. 
Throughout this paper we use the small scale formalism, so that instead of spherical expansion we work with plane 
wave expansion. This simplifies the expressions and reduces the computational time of simulations. The generalization 
to all sky coverage is presented in . 

The observed CMB temperature in the direction 9 is T(0) and equals the (unobservable) temperature at the last 
scattering surface in a different direction, T(9 + 69), where 69 is the angular excursion of the photon as it propagates 
from the last scattering surface to us. In terms of Fourier components we have 

T(9) = f{9 + 69) 

= (2vr)- 2 J d 2 l e d < e+se} f{l). (1) 

To extract the information on the deflection field 69 we consider derivatives of the CMB temperature. If the CMB 
is isotropic and homogeneous Gaussian random field then different partial derivatives are statistically equivalent and 
their spatial properties are independent of position. Lensing will distort these two properties of the derivatives. The 
derivatives of the temperature field are to lowest order 

T a (9)^^-(9 + 69) 

= (6 ab + $ ab )f b (9 + 69), (2) 

where & a b = is the shear tensor and a,b = x, y. The components of the shear tensor can be written in terms of 
the projected mass density K and the shear fields 71 and 72, 

<$>xx + $yy = -2k 

2<S> xy = -2 72 . (3) 

Convergence and shear are related to each other through the Fourier space relations, 

7i(I) = k(Qcos(2&) 7 2 (J) = K (/)sin(2^) ! (4) 

where <pi is the azimuthal angle of the Fourier mode I. The convergence k can be related simply to a radial projection 
of density perturbation ( Q ) 

« = ^a» / g(x',x)-dx'- (5) 

2 Jo a 

Here \ is the comoving radial coordinate of last-scattering surface and r(\) is the corresponding comoving angular 
diameter distance, defined as K^ 1 / 2 sin tO-^x, X> ( — K)~ x / 2 sinh(— K) x ^ 2 \ for K > 0, K = 0, K < 0, respectively, 
where K is the curvature, which can be expressed using the present density parameter f2o and the present Hubble 
parameter Hq as K = (£Iq — 1)H 2 . The density parameter fig can have contributions from matter density fl m , as 
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well as from other components such as the vacuum density J7a- The radial window over the density perturbations 8 
is g/a, where g(x' , x) — r ^ x ^r(x)~ X ^ 1S a sna P e( i curve symmetric around x/2 and vanishing at and x, while a 
is the expansion factor. Note that for a flat Q m = 1 universe 8 oc a in linear theory and the weighting is symmetric 
around x/2, so that the peak contribution is coming from z = 3. 

The angular power spectrum of convergence is defined as {k(1)*k(1')) = 0**5^ and has ensemble average pi 

cr-uatHtf Jlpg)ft(*= 4*) *• (•) 

Here Pg(k,r) is the 3-d dark matter power spectrum that is integrated over the past light cone and is in general a 
function of time r and wavevector k. This equation has been derived using Limber's equation, which is only valid in 
the small scale limit. The general solution in terms of the line of sight integral over the spherical Bessel functions is 
given in |Q. 

We consider next the quadratic combinations of the derivatives in equation (j^) and express them in terms of the 
unlcnsed field to lowest order in the shear tensor, 

S=[T 2 X +T*] (9) 

= (1 + $o + $yy)S + {® xx - %y)Q + 2® xy U 

Q ee [T 2 X - T 2 } (9) 

= (1 + + ®yy)Q + ~ $yy)S 

U = 2 [T x T y ] {9) 

= (1 + ^ + $yy)U + <Z> xy S, (7) 

where S, Q,U are the corresponding quantities in the unlensed CMB field at 9 + 59. The notation used here makes 
the analogy with CMB polarization: Q ± ihi in equation (^) have spin ±2 just like the Stokes parameters used to 
describe the CMB polarization, while S is a spin quantity (a scalar) and is rotationally invariant. 

Equation (Q) shows that the measured S, Q and U are products of the projected mass density and shear with 
derivatives of the unlensed CMB field. Thus the power spectrum of 5, Q and U will be a convolution of the power in 
the CMB and that of the projected mass density. The general expression of this convolution is quite involved, so we 
will discuss it in the two limits where the expressions simplify considerably, the limits of large and small scales relative 
to the CMB correlation length £. The large scale limit is sufficient to analyze the potential of MAP and Planck future 
missions. For experiments with higher angular resolution the full convolution will be necessary. 

Throughout the paper we compare the result of our analytical estimates with those of numerical simulations. To 
simulate the lensing effect on the CMB we first generate on a fixed square grid a projected density map. We then 
Fourier transform the convergence and compute the displacement vector 89 by using the Fourier relation 

89 = 2i^n (8) 

and then transforming it back to real space. We next generate on a fixed grid of the same size as above a random 
realization of temperature field T, using an input CMB power spectrum. We construct the lensed temperature map 
do that for each point in the lensed temperature map we use the corresponding displacement vector to determine 
from what position on the original grid the photons came from. This position does not generally coincide with a grid 
point in the unlensed map so we use cloud-in-cell interpolation to compute the value of T in the original map at this 
position. Cloud in cell interpolation smoothes the field so relatively small grid sizes are needed to avoid this unwanted 
effect. We achieve this by increasing the size of the array when performing this step. 



A. Large scale limit 

We start with the reconstruction in the limit of large scales. In this limit we average over many CMB patches 
to detect the weak lensing signal. We begin by noting that in the absence of lensing the isotropy of the unlensed 
background implies that 

(S)cmb = crs 

(Q)CMB - 

(U)cmb - 0, (9) 
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where a s = (T£)cmb + (Ty)cMB = 2(T%)cmb = 2(T^) CMB . This average can be computed in terms of the CMB 
power spectrum Cf T , 



°s=yf^ci>. do) 

If the mean of the components of the shear tensor vanish we have (S) = (S). There are residual quadratic terms in 
the shear tensor that contribute to equation (|l^) , but they are negligible in most cases of interest because shear is 
expected to be small. 

The average of equation (0) over an ensemble of CMB fluctuations gives, 

{S) C mb = (1 - 2n)a 2 s 

(Q)CMB — — 271CT5 

(U)cMB = -272CT|. (11) 

The physical interpretation of these equations is simple: k will stretch the image, which makes the derivatives smaller. 
Its effect is isotropic and only changes the value of S. The shear produces an anisotropy in the derivatives, in the 
same way as it creates an ellipticity in the shape of a circular background galaxy. This can be extracted from the 
particular combination of the derivatives used here. 

We can reconstruct k by studying the statistics of S and by combining the shear obtained from Q and IA. It is 
convenient to change the variables to 

1 



so that, 
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s' = - 


h 






Q' = - 


Q_ 
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U' = - 






as 


{S')cMB 


= 2k 


(Q')cMB 


= 2 7 i 


(M')cmb 


= 2 72 



(12) 



(13) 

We will drop the primes in what follows, as we will only use these quantities for the rest of the paper. Note that it is 
necessary to have a quadratic combination of temperature field for an unambiguous reconstruction of the convergence. 
This means that any reconstruction will have noise arising from intrinsic fluctuations in the CMB even in the absence 
of detector noise. Below we will quantify this intrinsic CMB noise. 

In the following we use the formalism developed to characterize the CMB polarization field pj] . We can combine 
Q and U to form £ and B, two spin zero quantities, which in Fourier space are defined as 

£(l) = Q(J)cos(2(^)+W(Z)sin(2<^) 

B(l) = Q{1) sin(20j) - U(l) cos(20 z ). (14) 

from which the real space £{0) and B(9) can be obtained by Fourier transformation. Equivalently these can be 
constructed directly from quantities in real space 

£(0)= J d 2 9' w(|0'-0|) Q r (d') 

B{6) = J d 2 6' uj{\0' - 0\) Ur(O'). (15) 

We have defined Q r and U r , the derivative shear in the polar coordinate system centered at 9. If = then 
Q r = cos 2(f)' Q{0') - sin 20' U{0') and U r = cos 20' U{9') + sm2(f>' Q{9'). The window is w(0) = -1/tt6» 2 (0 ^ 0), 
u)(0) = (0 = Q). From equations (|l|) and ^ it follows that 
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(S)cmb = 2k 

(£)CMB = 2k 

(B)cmb = 0. (16) 

Equation ( [lif ) shows that the average of 5 and £ can be used to reconstruct the projected mass density. However, 
in any particular direction on the sky the CMB derivatives can take any value so there is an intrinsic noise in the 
reconstruction coming from the random nature of the CMB. 

To describe the CMB noise we need to calculate correlation functions between S, Q and U. These can be expressed 
in terms of the correlation functions of the derivatives of the unlensed CMB field. For simplicity we consider two 
directions separated by an angle 9 in the x direction, 



c xx (6) = (f x {o)f x (e)) C MB 

= (2^)- 2 J d 2 l e «-««»*i l 2 Cf f cos 2 , 

= j l ^i 2 cr[Mie)-Mie)] 



Cyy(9) = (fy(Q)f y (9)) CMB 

= (2tt)- 2 J d 2 l e «-««»*i pQ ff sin 2 ^ 

= J^^r [Mm + Mm 



l[Co(0) + C 2 (9)} 



C xy (9) ee (f x (0)f y {9)) CM B 

= (2tt)- 2 J d 2 l e il - ecos <t» Pep cos ^ sin ^ 

-0, (17) 
where Cq(9) and C2(6) are defined as the integrals over l 3 Cidl/2ir weighted with Jq(W) and ^(W), respectively. The 



real space correlations of S, Q, U are 








N ss (9) EE 


(5(0)5((9)> 


= HC 2 XX + C 2 yy )/a 2 = 


{Gl + Ct)/a% 


N QQ (9) EE 


(Q(o)eo?)) 


— 2{C XX + C yy )/ a% = 


(C 2 + C 2 )/a 2 


N uu (9) EE 


(U(0)U(9)) 


— iC xx C yy / '<j s = (C 2 


-Cl)/a 2 


N SQ (9) EE 


(S(0)Q(9)) 


= ^^PxX ~ Cyy)/CT S = 


2CoC 2 /os 


N su (9) EE 


(S(0)U(9)) 


= 




N QU (9) EE 


(Q(0)U(9)) 


= 0. 





(18) 

With the normalization we have chosen the correlations for S, Q and U at zero lag are equal to one. The correlations 

18|), in a similar way 



when the separation is not along the x axis can be obtained by rotations of those in equation ( 
as done for the correlations of the Stokes parameters that describe the CMB polarization |l^ , |l"3 [ 

To reconstruct the power spectrum of k it will be necessary to have expressions for the noise power spectra Nf , 
Nf £ and Nf £ , defined as 

(W(l)W(l)) = 4Q KK + 7V ; WVV ', (19) 

where W stands for S or £ . For B only the noise term iVf contributes. These power spectra can be obtained from 
the real space correlation functions using Q 

Nf s = 2tt / 9d9 N ss (9) J o (Z0) 
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ATf £ =ir J 6d9{[N QQ (d) + N uu {6)] J a {l9) + [N QQ (d) - N uu (6)} J 4 (W)} 
iVf B =tt J BdS {[N QQ (9) + N UU (9)] J (16) - [N QQ {9) - N uu {6)] J 4 (W)} 
Nf £ = 2tt J 9d9 N S£ {9) J 2 {W). (20) 
From equations (Q~8h and (j20|) it follows that the CMB noise power spectra are, 



N ss 



2ir f 

— / MQ{Cl + Cf)Jo(W) 



jy£fi = / 0d6[C$J o (W) + ClJi(W)] 
a s 

N? B = % \ edd[ClJ (W) - Cf J 4 (Z0)] 



S 



ivf £ = / 9d9C C 2 J 2 {W). (21) 

The top panel of figure [l] shows the correlation functions from equation (18), while the bottom panel of the same 
figure shows the different noise power spectra using a standard CDM temperature power spectrum. The correlations 
in the CMB derivatives drop significantly for angles larger than £ = 0.15°. The power spectra shown in figure ^ 
demonstrate that the large scale behavior of the correlations is like white noise. In the limit of low I we can obtain 
the power spectrum from equation ( |2l| ) by integrating over angle first and using the orthonormality relation of the 
Bessel functions, 

limiVf s = 2tt J 9d9{C$ + Cl)/a 2 s 

= (2nal)- 1 J 9d9 J l 3 dl J l' 3 dl'Cf 7 ' C$ T [J Q (16)J Q {1'6) + J 2 (19)J 2 {1'9)] 

I fl 5 dl(C?T) 2 



.N? s 



7T (Jl 2 dlC? T ) 2 

\imN? £ = 2n [ 6ddC$/al = lim Np B = -lhn. 

i^o 1 J 0/ s i^o 1 2 i^o 

lim Nf £ = 0. (22) 

We can use these to define the correlation length more precisely as £ 2 = Nf s which gives £ ~ 0.15° used above. It is 
interesting to note that although S, Q, hi all have the same variance the power spectrum of S is twice that of £ or B. 
£ and B have a different correlation length from <S, which reflects the spin nature of the shear variables. This leads 
to £ giving a factor of 2 higher signal to noise than S in the reconstruction. 

In order to get an accurate measurement of the projected mass density or the shear we need to average over several 
"independent" patches of size £ 2 , each of which has a variance of order unity. This sets the basic requirement that 
has to be satisfied if the reconstruction is to give a positive signature. The signal to noise in each patch of size £ 2 is 

— 1/2 

S/N ~ 2k. This means that we need to average over N tch ~ 2\k\ to get a signal to noise of order one. The issue 
is whether on scales larger than £ the weak lensing signal is sufficiently strong to be detectable. Even though CMB 
is sensitive to matter density fluctuations up to z ~ 1100 and the RMS convergence is significantly higher than for 
galaxies at say z ~ 1, it is still well below the CMB noise if £ = 0.15°. Only with sufficient small scale power can 
individual structures be reconstructed with high enough signal to noise. This is discussed further in the following 
sections. 

Equations ( pl| ) are needed to asses the noise for the reconstruction of K and to subtract the noise contribution to 



its reconstructed power spectrum in equation (19). To compute the variance of the power spectrum reconstruction 
we also need to know the RMS of the CMB noise power spectrum. The CMB noise is a fourth order statistics of 
CMB field. Fortunately in the large scale limit it can be considered Gaussian. This is a consequence of the central 
limit theorem: the long wavelength modes of S, £ and B are obtained by adding a lot of independent patches making 
them effectively Gaussian. 
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FIG. 1. The upper panel shows the correlation function of S, Q, U for SCDM. The lower panel shows the power spectra of 
SS, ££, BB and S£. 



If the CMB noise can be considered Gaussian then the covariance matrix Cov{XX') for the noise power spectra 
can be expressed in terms of these [[ll] , 



Cov[(Af s ) 2 ] = 
Cov[(iVP) 2 ] = 



Cov[(iVf £ ) 2 ] 



Cov{N? s Nf s ) 
Qov{N? s N? £ ) 
Cov(Nf £ N? £ ) 



21 + 
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2 




21 + 
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1 




21 + 
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2 




21 + 
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2 




21 + 
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2 




21 + 


1 



-Kn^Y + n^n; 



N S£ N ss 



N? £ Nf £ . 



(23) 



In the large scale limit Nf £ < Nf £ and Nf s 



2Nf £ 



2Nf B . 



The covariance matrix becomes diagonal with 



Cov[(N l ££ ) 2 ] = Cov[(iVf 8 ) 2 ] = Cov[(Nf £ ) 2 ] = Cov[(7Vf s ) 2 ]/4. In figure |the covariances obtained in the simulations 
are compared to those derived under the Gaussian approximation. The agreement is excellent for I < 1000 and only 
seriously breaks down beyond that for N S£ . We have also verified with Monte Carlo simulations that the off-diagonal 
terms are negligible compared to the diagonal ones. 

In equation ( |l9| ) we have three estimators for the convergence power spectrum, SS, SB and the cross correlation 
S£. Using the covariance matrix in (|23| ) we can construct a minimum variance combination of the three estimators. 
In the large scale limit the covariance matrix is diagonal and we can weight each of the estimators by the inverse of 
noise variance, 



:(Ci 



ss 



n; 



SS\ 



Nf £ ) + -(Cf £ 



Nf £ ) 



(24) 
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FIG. 2. Simulation results for the covariances of N t , N t N t and N t normalized to their values for Gaussian noise 
given in equation (p3|). The Gaussian approximation is an excellent one for I < 1000. 

Note that we have assumed in equation (|23|) that only the CMB noise contributes to the variance of C ; KK estimates. 
There is also the cosmic variance contribution 2nr[(G/ reK ) 2 but we will be show in the following sections that this term 
is significantly smaller than the CMB noise term and can be ignored, unless the CMB has more small scale power 
than expected within the currently popular cosmological models. 



B. Small Scale Limit 



So far we have described the gravitational lensing effect on large scales compared to the correlation length of 
the CMB. In this limit we average over (almost) independent patches of CMB fluctuations, which is similar to the 
procedure used in weak lensing, where we average over independent ellipticities of background galaxies. We now turn 
to scales much smaller than the correlation length £. In this limit we can take the derivatives T x and T y to be constant 
across the field. The physical picture is very different from the one discussed in previous subsection. Weak lensing 
in this limit acts as a generator of small scale power. To see this we can imagine the effect of a small clump of mass 
on a pure gradient of temperature field. A mass clump will magnify and stretch a small patch, which will change 
the slope of the gradient at the position of the clump. The resulting temperature field is no longer a pure gradient, 
but has a small wiggle superimposed on it. Small scale power has been generated only where the gradient of primary 
anisotropics is non-zero because where the surface brightness is constant its conservation requires that it remains so 
even in the presence of lensing. 

We now attempt to use this physical picture to reconstruct the weak lensing signal. Introducing a = T x jas and 
b = T y /as we have 



S = (2k - l)(a 2 + b 2 ) + 2 7i (a 2 - b 2 ) + 2 72 lab 
Q=(2k- l)(a 2 - b 2 ) + 2 7l (a 2 + b 2 ) 
U = (2k- l)2a& + 2 72 (a 2 + b 2 ). 



(25) 



We did not subtract unity out of the definition of S as in equation (O), because even though on average ((a 2 
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b 2 ))ci\iB = 1, we are considering a very small field over which a and b are approximately constant with (a 2 + b 2 ) ^ 1 
in general. 

We want to determine a, b, k, 71 and 72. If we assume that the mean value of k, 71 and 72 over the field is zero 
then 

(5) = a 2 + 6 2 
(Q) = a 2 - b 2 

{U) = 2ab, (26) 

where the mean value is taken over a region large enough that the mean shear and projected mass density vanish but 
over which the CMB gradient remains nearly constant. In practice this can be achieved by filtering out the small 
scale power, so that the remaining power is largely dominated by primary anisotropies. We can then compute a and 
b and use them in equation (^6|). 

Once a and b are determined it would appear that we can determine the convergence and the two shear components 
by solving 

5 \ / (5) (Q) (U) 

C = (Q) (S) 1( -2 7 i 1 . (27) 
U ) \ (0) (U) (S) 

However, the determinant of the matrix above vanishes so we cannot determine convergence and shear independently. 
Instead we have to use the relations between shear and convergence to reconstruct the signal. In Fourier space, for I 
large enough that we can consider a and b in equation (^5|) as constants. We get, 

S{1) = 2k(1) [(a 2 + b 2 ) + (a 2 - b 2 ) cos(2<^) + 2o6sin(2^)] 
Q(l) = 2/c(I) [(a 2 - b 2 ) + (a 2 + b 2 ) cos(2^)] 

U{1) = 2k(1) [2ab + (a 2 + b 2 ) sin(2^)] . (28) 




We can compute £ and B, 

8(1) = Q(i)cos(2^) +U(l)aia(2<t>i) 

= 2k(1) [(a 2 + b 2 ) + (a 2 - b 2 ) cos(20 ; ) + 2a6sin(20 z )] 
= S(l) 

B{1) =W(i)cos(2^) - Q(i)sin(2^) 

= 2k(1) [2abcos(2^) + (a 2 - 6 2 ) sin(20 z )] . (29) 

We see that the estimators of the convergence are not independent, £ (I) = S(l), as argued above. 

To obtain an estimator of convergence we can use any of the expressions in equation ( pq ) . This requires dividing 
with the combination of derivatives of T and will be very noisy if both a and b are close to 0. This is a consequence 
of the fact that lensing cannot generate power where there is no gradient. Such a reconstruction will therefore have 
variable noise. One solution to this problem is to filter the map with a variable filtering length. Note that we have 
assumed a and b are constant across the map, while in reality they will change as well. This means that in practice we 
have to divide the map into chunks over which the long wavelength modes are approximately constant, or use more 
sophisticated methods such as the wavelet analysis. 

The procedure outlined above is quite involved and we did not attempt to implement it in this paper. It simplifies 
considerably if we are only interested in the power spectrum. In this case we may still use S(l) and B(l) to get, 

Cf s = 4Cr [(a 2 + b 2 ) 2 + {a 2 - b 2 ) 2 cos 2 (2^) + (2a6) 2 sin 2 (2^) 

+ 2(a 2 + b 2 )(a 2 - b 2 ) cos(20j) + 4(a 2 + 6 2 )a6sin(2^) + 4(a 2 - 6 2 )a6cos(20j) sin(20j)] 

Cf 8 = 4Cf K [+(2afe) 2 cos 2 (20j) + (a 2 - 6 2 ) 2 sin 2 (2^) 

- 4(a 2 - 6 2 )afecos(2^) sin(20j)] . (30) 

If we average over ^ and use ((a 2 +b 2 ) 2 ) c M b = (f^+f^+2f 2 f 2 ) C MB/<J 2 s = 2 and ({a 2 -b 2 ) 2 ) CMB = ((2ab) 2 ) CMB = 
1 we find, 
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Cf s = 12C ; KK 

Cf B = 4C*f K . (31) 

Note that there is no need to use local estimates of a and b to get an estimator of the power spectrum, as long we have 
a CMB map with enough uncorrelated patches so that the quadratic combinations of a and b can be replaced by their 
averages. In this limit the power spectra of <S and B again directly give estimates of convergence power spectrum, 
without any convolution and also without noise contribution from intrinsic CMB anisotropics (although there will be 
contribution from instrumental noise). This is because in this limit all the power is generated by gravitational lensing. 
It is important to note that in this limit B does not vanish, but actually gives an estimate of the k power spectrum, 
while power in S is a factor of 3 bigger than power in B. We will show below that these predictions are well recovered 
in the Monte Carlo simulations. 



III. RECONSTRUCTING THE MASS PROFILE OF A CLUSTER 

To illustrate our method we first apply it to reconstruct the mass profile of a massive cluster of galaxies. In order 
for the large scale limit method to work we need the CMB to vary on scales smaller than the cluster, so that the 
effect of lensing can be measured after averaging over independent patches. Since we do not expect the primary 
CMB anisotropies to have sufficient power on arcminute scales we use the Ostriker-Vishniac (OV) effect instead in 
our example using the power spectrum from (lj] as an estimate. This is not necessary the only possible source of 
anisotropies on small scales and there may be other sources as well, including the primeval galaxies and QSO emitting 
in IR and radio. As long as there are many sources in the beam and their redshift is higher than the cluster redshift 
the analysis is identical to the one presented here. The only difference is that these sources are not at z ~ 1100 but 
at lower redshift so the window function g in equation (ph changes relative to the one for z ~ 1100 sources. If the 
redshift distribution of the sources is not known then we may in principle use lensing on a cluster with known mass 
profile to determine it. For individual sources that are resolved it may be better to use their shapes as an estimate of 
the local shear, which is the usual procedure in the case of galaxy lensing. 

For the method to be useful in practice the small scale anisotropies have to be above the detector noise. This could 
be achievable with future interferometric or bolometer experiments with long integration times on a small area of the 
sky, reaching fj,K sensitivities on subarcminute pixels. Another complication in this case relative to our analysis in §2 
is that the distribution of anisotropies is likely not to be Gaussian. This complicates the issue of noise analysis and 
statistical significance of the results, but does not affect the average of the reconstruction. We will ignore these issues 
in the following, since we are mainly interested in having a simple example with which we can test our method. 

Figure Sa shows the simulated input cluster, with the profile chosen as a projection of a truncated isothermal sphere. 
In figure |[b we show an unlensed simulated CMB field, assuming OV type of power spectrum. In figure ||c we show 
the lensed CMB map. We see that the cluster has magnified the CMB in the center, so in the absence of small scale 
power the central region is smoother. In figure |§|d we show the shear as measured by the CMB derivatives (Q and IX) 
and in the background the convergence field S. 

Both S and £ are estimators of the convergence, but figure ||c illustrates a possible advantage of the reconstruction 
based on shear over that based on S in the case of clusters. The shear based reconstruction is nonlocal as shown from 
equation (^). The shear from the whole map is used to reconstruct the density at the center. Thus we are averaging 
over many different independent CMB regions. On the other hand the reconstruction based on S is local. It will not 
work very well in the center because the magnification induced by the lensing will expand a small region so unless 
there are additional small scale perturbations which are magnified there will be fewer independent regions to average 
over. This effect can be seen in figure ||c. This is of course an issue only for clusters that are close to being critical 
and the magnification is no longer small compared to unity. For other, more linear structures the two reconstructions 
give similar results, although S is still noisier than £ by a factor of 2 as derived in equation (B2|) . 

To be more quantitative we show in figure ^ the reconstructed projected radial mass profile from the shear and the 
convergence. We also show the result of the reconstruction procedure on the unlensed CMB field to illustrate the 
noise level, which for this case with a lot of small scale structure is negligible. The reconstructed profile tracks well the 
input profile outside the center, but falls below it in the center because of insufficient number of independent patches 
there. In our reconstruction algorithm we are forcing the mean k to be zero and because the reconstruction of the 
center of the cluster is systematically lower due to the averaging, the reconstruction is also systematically lower on 
the outside of the cluster. As expected this is more pronounced for S than for £ reconstruction because the former is 
noisier and relies more on poor quality information from the center of the cluster. These effects can be calibrated from 
the simulations and overall the method can reconstruct the true projected mass density outside the core, assuming 
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FIG. 3. Top left panel: input cluster on a 6'x6' field. Top right panel: unlensed CMB map. We assumed that the 
Ostriker-Vishniac effect could be detected with sufficient signal to noise to be used in the reconstruction of ft. Bottom left: 
CMB field after being lensed by the cluster. Bottom right: the background shows the 5 field while the rods represent the shear 
variables Q and U, both of which can be used to reconstruct the density profile. 
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(arcseconds) 

FIG. 4. Reconstructed radial profile using the data in Figure 2. The points marked with £ (S) correspond to the recon- 
struction of k based on £ (S). The points marked unlensed correspond to the result of applying the £ reconstruction to the 
unlensed field. The input k was taken to have zero mean. 

that the necessary conditions discussed above (measurable primary anisotropies on scales smaller than the typical 
cluster scale) are satisfied. 

IV. RECONSTRUCTING THE POWER SPECTRUM OF PROJECTED DENSITY 

A. CMB Noise 

In this section we investigate the power spectrum reconstruction of convergence using the method developed in 
§2. The power spectrum can be reconstructed from £, S or their cross correlation. We will first discuss the power 
spectrum reconstruction in the absence of detector noise and assuming infinite angular resolution. In this simplified 
case the only source of noise is the CMB itself. We do not add small scale secondary anisotropies in this section, 
so the results should be representative for primary fluctuations and any additional small scale fluctuations at high 
redshift would further improve the performance. In figure ||a we show the power spectra for the unlensed CMB field 
(the CMB "noise") calculated using equation ( pl| ) and the result of our Monte Carlo simulations. The agreement 
between analytic prediction and simulations is very good. We recover the white noise behavior of Nf s ,Nf £ and N® B 
on large scales and also the relatively small amplitude of the cross correlation Nf £ . 

On large scales these power spectra are to be compared to 4C ; KK (figure ||a) . For the purpose of illustrating our 
method in this paper we adopted the convergence power spectra of the "concordance" model JI^] fl m — 0.3, £1a = 0.7, 
T = Q m h = 0.2, n = 1 and erg = 1. It can be seen from figure ||a that on large scales the power of the convergence 
is small compared to that of the CMB noise. As a result we can only hope for a statistical detection of the lensing 
effect and the method will not allow to make a map of k directly unless there is more small scale power in the CMB, 
as discussed in §3 in the context of cluster reconstruction. Even if k reconstruction has a low signal to noise its 
power spectrum may nevertheless be recovered with high statistical significance, because we can average over many 
independent patches on the sky. This is seen from the difference between the power spectra in the lensed and unlensed 
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FIG. 5. (a) Nf s , Nf £ , Ni and Nf £ power spectra for the unlensed CMB field. The lines correspond to the results of 
equation (pit). Also shown is the weighted average of the recovered power spectra from S, £ and their cross correlation together 
with the input 4Cf K . (b) Ratio of the recovered power spectra to the input Cf K . The signal in B is also shown. 
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field. The weighted average of the three estimators (equation gj) is also shown in figure [|a together with the input 
convergence spectrum. The intrinsic CMB noise is about 10-20 times higher, yet the difference agrees well with the 
input k power spectrum, showing that there is a significant statistical detection of the signal. 

The reconstruction of the power spectrum is further studied in figure ^d, where we show the CMB noise subtracted 
power spectra of 55, EE and BB and the cross correlation spectrum SE divided with the input matter power spectrum. 
On large scales the reconstructed spectra of SS, EE and SE all give unbiased estimators of 4C ; KK so the average of the 
output to input ratio is 1, while that of BB if consistent with the pure noise and the ratio averages to 0. The latter, 
although not giving an estimate of the signal, can provide an important consistency check on whether the signal seen 
from £ and S is indeed of cosmological origin or not. It may also help identify the cosmological part of the signal if 
additional contributions are present, such as in the case of non Gaussian CMB fluctuations or contamination from 
foregrounds. 

The large scale limit is valid up to I ~ 1000. On very small scales, roughly I ~ 5000 and beyond, the reconstructed 
power spectra from £ and S give 12C ; KK , while that obtained from B gives ACf K so that the ratios are 3 and 1, again 
in agreement with analytic prediction in equation ([}l]). At very high I CMB noise is negligible, because there is no 
power present in the CMB itself on those scales. 

In the intermediate regime between I = 1000 — 5000 neither small scale nor large scale limits apply and the 
reconstruction becomes a complicated convolution of CMB anisotropies and lensing signal. Information on the power 
spectrum can still be obtained even from this regime, but it is better to approach the reconstruction in a parametric 
form, by parameterizing the power spectrum with a few free parameters that can be estimated by fitting the simulated 
spectra to the data, rather than by direct inversion. 

Because CMB noise is larger than the signal on large scales we may worry that a small error in the noise estimate 
would lead to a large error in the estimated power spectrum. This is most worrisome for the SS reconstruction and 
least important for the cross correlation SE because Nf s = 2Nf £ 3> Nf . The latter is even smaller than 4C ; KK 
on large scales so we can obtain an accurate estimate even without any noise subtraction. Fortunately the large 
scale behavior of Nf s , Nf £ Nf B and Nf £ is quite insensitive to the details of the power temperature spectra. For 
example we may worry that the CMB power spectra we measure will be affected by lensing or that we can only 
know the lensed CMB temperature spectra with some minimum scatter, limited by cosmic variance. On large scales, 
these effects change the mean of Nf t by approximately 0.5 % and add a 1% scatter to it. It will depend on the 
noise amplitude whether this is a significant source of error. For MAP these errors causes a shift in the amplitude 
comparable to the input power spectrum, while for Planck it is significantly below it and can be ignored. For MAP 
this problem can be solved by using B reconstruction. Because the effect on large scales is mostly an amplitude shift 
we may use the amplitude of Cf to make it consistent with 0, which would also properly determine the amplitude 
of the other spectra. 



Next we consider the influence of finite angular resolution and detector noise. In order to measure the convergence 
power spectra using the small scale limit the experiment must have enough angular resolution to probe these small 
scales. This effect of finite angular resolution is straightforward to include and we will not discuss it further, because 
future CMB satellites such as MAP and Planck will not be able to probe this limit. Here we discuss the more 
important effect of finite angular resolution and detector noise on the large scale limit reconstruction. It is easy to 
see that finite angular resolution has an effect on the intrinsic CMB noise amplitude. Because of finite resolution and 
detector noise small scale CMB power cannot be resolved and this leads to a larger correlation length £. This means 
that we have fewer independent patches to average over the CMB noise and the overall level of noise power spectrum 
is higher. Finite angular resolution also has a more subtle effect by reducing the transfered power. To understand 
it quantitatively we first compute the derivatives of the temperature field in the presence of some filtering function 
F(9), which involves the experimental beam and any additional filter we may want to use in the analysis 



B. Finite Angular Resolution and Detector Noise 




/ 
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+ (2TT)- 4 J d 2 l J d 2 qF(\l + q\)T b {l)^ ab {q)e l( - l+( l> 6 '. (32) 

To obtain the last expression we expanded F, $ ab and T a into a Fourier series and integrated over angle 0' . The large 
scale effect of the filter function can be read from the average over CMB of the square of equation (p2|), 



{TaT b ) C MB = 5 -fJ l ^F\ircr 

+ (2tt)- 2 J d 2 q<P ab (q)e^ d (2^)~ 2 J d 2 lF{l)F{\l + g|)Q 



TT 



(33) 



In the absence of filtering the integrals involving Cf T in both terms of equation ( |33|) are identical. This allows us to 
reconstruct K from the second term, while the first term gives as as in equation (|ll|). To describe the effect of the 
beam smearing we introduce a window function as the ratio of the two integrals 

= fd 2 iF(i)F { \i +q \)i 2 cr (34) 

which is in general less than unity and in terms of which we have 

(S(l)) = (£(l)) = 2k(1)W(1). (35) 

The effect of filtering is important even for the reconstruction of large scale k modes. This can be understood 
by looking at equation (|32|). The shear tensor modulates the amplitude of the temperature derivatives, so the 
information about a particular q mode of k is encoded as sidebands of the different I modes of the temperature. We 
are recovering the information back at the corresponding q by squaring the field, which appropriately combines back 
all the sidebands from every I to the correct q. Finite angular resolution is important even for large scale modes 
because the information about these modes in encoded mainly in the sidebands of the small scale temperature modes, 
which are strongly affected by the filtering function. 

To minimize the effect of beam smoothing we may filter the temperature before squaring it. The smallest effect is 
achieved if the filtering function F(l) is a constant, in which case (in the absence of a cutoff in /) there would be no 
effect. For this reason we choose to filter the temperature with the inverse of the beam e H'+ 1 ) CT i, j where a 2 is the width 
of the Gaussian beam, thereby reversing the smoothing effect of the beam. We can only do this up to a maximum 
l cu t (the value of which is determined below) not to amplify the detector noise on small scales. Thus our effective 
filter F(l) is equal to one for I < l cut and zero after that. Figure ^a shows two examples of W(q) for different values 
of l C ut, which correspond to minimum of noise for MAP and Planck experiments. Note that regardless of the form of 
the filter equation (|mJ) implies that W(q) — > 1 as q — > 0. Other forms of filtering are possible as well, but we show 
next that this filter comes close to being optimal and is easy to implement in the analysis. 

The contribution of detector noise on large scales to the different power spectra can be obtained using equation 
@ with Gj T Cj T + N TT . We take the detector noise to be white noise with constant power spectrum N 
for a range < I < l m ax- The power spectrum N TT is related to the noise level per pixel, a 2 = N TT l^ax I ^ n with 
^inax ~ NpixAw/n and a is expressed in the same units as T (e.g. fJ,K). If we want the correct normalization for the 
recovered Cf K we need to divide by the correct as, which has to be computed only with the CMB power spectrum, 
without the detector noise term. With our filtering prescription the large scale amplitude of the Nf power spectrum 
in equation (|22]) becomes 



( J^ ut ldll 4 F(l) 4 (Cr + ^ l+1 ^N TT ) 2 

(jidu 2 F 2 (i)cr) 2 



To derive the optimal filter F(l) we minimize Nf £ /W 2 (l )Cf K with respect to F 2 (l). In the limit I — > W(l) = 1 and 
this becomes equivalent to minimizing Nf £ in equation (|3^). Taking derivatives with respect to F(l) 2 and setting to 
we find 

C TT 

F{1) = Z2( C TT + e /(/+l)<x? 7V TT)2' (3?) 

In the large scale limit this gives F(l) 2 — (PCf 7 )^ 1 , which is roughly a constant for spectra that are close to scale 
invariant. For large / the noise term dominates over the CMB signal term and the filtering function F(l) goes to 0. 
The transition occurs where 
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FIG. 6. The upper panel shows the window function W 2 (l) as a function of I for MAP noise and beam properties where 
the appropriate value is l cu t = 600 and for Planck with l cut = 2000 (bottom panel). For the latter the window has no effect 
at I < 1000. Bottom panel shows amplitude of Nf £ in the limit I — > as a function of l cu t for the specifications of the MAP 
and Planck satellites, justifying l cut = 600, 2000 for the two experiments. Also shown is the no noise infinite angular resolution 
curve (Ideal), which levels off at I ~ 2000, showing that Planck is close to optimal in the large scale limit, because it measures 
most of the power in the CMB. 



Cj T ~ e l ^<N TT . (38) 

To compare this filter to the simple constant filter in figure [|b we show Nf £ as a function of l cut for the constant 
filter and the specifications of noise and angular resolution of MAP and Planck, as well as in the absence of noise 
and beam smoothing. In the absence of detector noise larger l cut is always better because the correlation length is 
reduced, but it saturates beyond I ~ 2000 for this model because there is very little CMB power on small scales. 
Once detector noise is included it eventually dominates at high l cut and we are better off not including these modes 
as they are mostly noise and do not contribute to the signal. In between there is a minimum which determines lcut 
depending on the noise and beam of the experiment. The amplitude of noise for this filter can be compared to the 
one in equation ( |37| ) which minimizes the noise in the large scale limit. The two give very similar results in the large 
scale limit. Because constant F(l) minimizes the effect of the window at higher / it is likely to perform even better 
than the filter in equation (p7p, which was derived under assumption W(l) = 1. For this reason we choose to adopt 



the simple constant filter instead of the one in equation (37) 



C. Relation to the Density Field 



We have shown that the quantity we are able to extract from the distortions of the CMB is the power spectrum 
of the projected density field weighted with a window g/a (equation ||). Since the more fundamental quantity is the 
3-d density field described by the power spectrum we would like to know the relation between the two power spectra. 
This relation is shown in figure |^a, where a logarithmic contribution to a given I mode as a function of 3-d wavevector 
k is plotted for a representative sample of models. These windows are relatively broad functions of k and peak at 
A = 2n/k = lOOO/i^Mpc for I = 10 and A = 30/i _1 Mpc for I = 1000. The exact value depends on the model. We are 
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FIG. 7. Left: logarithmic contribution to Cf K as a function of k for / = 10, 100, 1000 (the normalization is arbitrary). The 
models are flat CDM model (dotted), open CDM model with fl m = 0.3 (dashed) and cosmological constant model with with 
fl m = 0.3 (solid). All the models have F = f2 m /i = 0.21. Right: logarithmic contribution to Cf K as a function of 1 + z for the 
same models as above. 

therefore probing the power spectrum over a range of scales which extends to scales larger than any other method 
that directly traces dark matter. 

The next question we want to address is the redshift distribution of the contribution to a given Cf K . Even though 
the g I a window peaks at z = 3 for a flat universe this does not mean that the dominant contribution comes from this 
redshift. The relation between / and k depends on the shape of the power spectrum (equation ^). For any given I 
there is a range in k that contributes (figure Qa). The contribution from a high k mode comes from structures which 
are relatively closer to the observer (and so at lower z) than the structures that dominate the contribution for a lower 
k mode (equation ^|) . On large scales the matter power spectrum has a turnover so low k modes have very little power 
which implies that their contribution is smaller than that of higher k modes from more nearby structures despite 
the geometrical factor that suppresses these higher k modes. Thus for low I the contribution will be dominated by 
low z structures. The logarithmic contribution to a given I mode as a function of 1 + z is shown in figure and 
confirms these expectations. At smaller scales the power is more evenly distributed as a function of k and the peak 
contribution moves to higher rcdshifts. For / ~ 1000 it peaks between z ~ 2 — 3 with a long tail extending to higher 
z. At these scales we are therefore directly probing dark matter clustering at high redshifts. 

The two characteristic features of lensing on the CMB are the large scales and early epochs we can probe. In 
combination these facts guarantee that the power spectrum will be dominated by linear contributions. Figure ^ 
shows the difference between the Cf K calculated using the linear and nonlinear matter power spectrum. The spectra 
are shown as a function of I for the same set of models used in figure [?]. The nonlinear matter power spectra were 
computed using the linear to nonlinear mapping Jl8| |. Up to I ~ 1000 the power spectrum is dominated by linear 
modes. This has the advantage of making the results simple to interpret in terms of cosmological models, since no 
nonlinear corrections are necessary. Since MAP and Planck reconstruction in the large scale limit will not extend 
beyond I ~ 1000 they will be dominated by linear scales. At higher I nonlinear corrections become important. While 
this has the disadvantage that the interpretation becomes more complicated it has the advantage in that the power is 
boosted by almost an order of magnitude compared to the linear scales and so the lensing effect becomes more easily 
observable. 

V. CROSS CORRELATION WITH OTHER MAPS 

We have shown that individual structures cannot be reconstructed with sufficient signal to noise, unless the CMB 
has more small scale power than expected from primary anisotropies. One way to obtain a positive detection discussed 
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FIG. 8. Power spectrum C* K as a function of I for the same models as in figure [?]. Thin lower curves show the linear power 
spectra, thick upper curves the nonlinear spectra. 

in previous section is by combining information from independent patches of the sky into a measurement of the power 
spectrum. In this section we discuss another way to enhance the signal by using the cross correlation with another 
map with higher signal to noise. A cross correlation between signal induced by lensing and some other tracer has been 
explored before |l9| . The difference to our method is that they do not reconstruct k directly, but compute the photon 
deflection angle instead for which the 3-d matter distribution is needed. It can only be obtained from a redshift 
survey (such as SDSS or 2dF), under the assumption that light traces mass. Our method can also be used when 
such a 3-d distribution is not available. For example, we may try to cross correlate the reconstructed convergence 
map with Sunyaev-Zeldovich (SZ) map, which traces the integrated pressure along the line of sight. This would give 
positive cross correlation because of clusters, which should contribute to both convergence and SZ effect. Thus even 
if clusters cannot be individually detected using the method described in §3, they may still be detected statistically 
through this cross correlation, which would give information about pressure to dark matter ratio in clusters. Other 
possibilities for cross correlation are with X-ray background which is believed to trace large scale structure to z ~ 5 
and with galaxy catalogs from SDSS and 2dF, both of which would give information on how light traces mass on 
large scales. Yet another possibility is the cross correlation with the CMB itself, which would give p ositive detection 
whenever there is a significant contribution from the time-dependent potential in the CMB M]20|]21|l . To test the 
most optimistic possibility we may estimate the signal to noise of the cross correlation when we have another perfect 
map of convergence. For any of the probes mentioned above the actual signal to noise will be lower, since they 
will only partially correlate with the projected mass density. It is nevertheless of interest to see whether such cross 
correlations would be useful at all given such a low signal to noise on individual structures in the maps of S and £. 
The cross correlations with the input k give 

Cf K = 0. (39) 

Figure |^ shows the cross correlations in equation ( |39| ) for a 100° x 100° field where the CMB is measured with 
negligible detector noise and infinite angular resolution (the results are almost the same for Planck specifications). The 
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FIG. 9. Cross correlation of the input re and the £ and S statistics. The results correspond to a 100° x 100° simulation 
observed with an ideal CMB experiment with negligible detector noise and infinite angular resolution. 

agreement is remarkable, proving that there is enough signal in the re maps recovered using our technique to get useful 
information about the convergence power spectra by cross correlating them with other maps. This is not surprising, 
since there are about 5xl0 5 independent patches of size £ 2 = (0.15 ) 2 and we need about (Nf £ /4C i ' tK ) 2 ~ 100 — 400 
of them to obtain signal to noise of unity. 

In reality most tracers of the underlying mass will not be perfect because they correlate only partially with re. In 
particular we assume that we have a map Y that correlates with re and has a cross correlation C ; Ky . Then the cross 
correlations between the W = S, £ maps and the Y map give, 

C ; wy = 2C? Y W{1) (40) 

where W(l) is the window defined in equation (pi[). The correlation of the Y map with B vanishes, Cf Y = 0. The 
covariance matrix for the two correlations is, 

Cov{Cf y C™ y ) = {Cj Y + N YY ){4C<rW 2 {l) + N™' w ) + ACfCfW 2 {l), (41) 

where the first term is usually dominant. Both the decorrelation and the noise reduce the overall signal to noise of 
the cross correlation. For example, in the case of cross correlation with the CMB only the largest modes are strongly 
correlated because the time dependent gravitational potential does not contribute to the CMB anisotropics on small 
scales. The overall signal to noise is significantly reduced compared to the idealized example discussed here, but is 
nevertheless detectable with future satellite missions for reasonable low density models [Jl0| . 

VI. CONCLUSIONS 

We have developed a method to reconstruct projected mass density from observed cosmic microwave background 
maps. The method consists of taking derivatives of temperature field and averaging their products. Particular 
combinations on average correspond to shear and convergence of gravitational lensing, which can be expressed as a 
line of sight integral over the density field. 
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We have presented three possible applications of the method. First we applied the method to simulated clusters 
showing that it can be successfully used to reconstruct their projected mass distribution if there is sufficient small 
scale power. For this to be successful we require small scale power beyond the one given by primary anisotropies, 
which could be provided by secondary processes such as Ostriker-Vishniac effect or primeval galaxies. The expected 
signal from larger structures, such as filaments and superclusters, is smaller and therefore more difficult to directly 
observe. 

A second and more promising application is to average the reconstructed map to extract the power spectrum of 
fluctuations. In simulations the method successfully reconstructs the input power spectrum both in the large scale 
limit (where we are averaging over many independent patches of CMB to reduce the noise) and in the small scale 
limit (where lensing generates additional small scale power). Future satellite missions should be able to measure this 
signal at least on large scales p5| and provide additional constraints on the cosmological parameters. Interferometers 
may be able to determine the power spectrum on smaller scales as well. This method is complementary and in many 
cases more robust than the more traditional methods of determining the power spectrum. In comparison to the power 
spectrum from galaxy surveys the main advantage is that there is no assumption on how light traces mass, which is 
still poorly understood at present. In comparison with the weak lensing surveys planned for the future the advantage 
is that the redshift distribution of the source is well known (z ~ 1100) and much higher than for any galaxy survey. 
It also does not suffer from intrinsic correlations between the galaxies, which may mimic the weak lensing signal. In 
comparison with other direct tracers of dark matter such as velocity flows and Ly-a forest the method presented here 
recovers the power spectrum over a larger range in scale and is less sensitive to basic assumptions underlying the 
method such as the error distribution of galaxy distances or assumptions of the IGM at high redshift. In addition, the 
method presented here also gives power spectrum information on much larger scales than reachable by other methods. 
These are likely to be linear and may allow us to deconvolve the power spectrum to obtain the 3-d power spectrum 
using the methods developed in p2fl . 

A third promising application of the reconstruction is its cross correlation with other maps that trace large scale 
structure. Some of these are X-ray background, galaxy surveys and CMB itself (both the primary anisotropies and 
SZ contribution). This again allows us to average over many independent patches to reach a positive detection and 
may give even higher statistical significance than the power spectrum if the two maps are well correlated. In this case 
we can learn not only about the dark matter clustering, but also about its relation to X-ray and optical light or about 
a time dependent gravitational potential. Overall, CMB fluctuations may hide a whole new information treasure in 
its pattern and its extraction would provide important information on the universe we live in. 
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